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1 INTRODUCTION 

Quantum Chromodynamics (QCD) [p], is regarded as the best, if not the only, viable 
theory of the strong interactions. Recent theoretical and experimental developments have 
significantly increased our ability to perform quantitative tests || and have deepened 
our understanding of hadronic interactions. This is particularly true for higher energy 
processes where the decreasing value of the coupling constant a s (fi) allows reliable results 
from a perturbative expansion ||-[]n|]. In addition, the non-perturbative transition from 
the fundamental objects in the theory (quarks and gluons) to observed particles (hadrons) 
has a smaller influence on measured quantities. At the energies now accessible, one is 
expected to be far away from the long distance regime where perturbation theory breaks 
down, and from the ultra-short distance regime where one might witness the onset of new 
dynamics. This, however, does not excuse one from vigilance for significant deviations 
from perturbative predictions. 

Tests of QCD in hadron-hadron collisions display a parallel development in both theory 
and experiment. The earliest measurements of high p t hadron production at the ISR 
belied the then hidden partonic component of the proton. The observation of jets and 
early tests of QCD at the SppS jn| were largely qualitative, yet they demonstrated the 
predictive power of the theory at leading order in perturbation theory. Currently we are 
entering a period where the emphasis is being placed on measurement precision. From an 
experimental standpoint this implies making measurements with high statistics and small 
systematic uncertainties. From a theoretical standpoint, this means calculating quantities 
at successively higher orders in perturbation theory, and using constraints from a number 
of sources (c.f. parton distribution functions) to pin down predictions. In the interplay 
between theory and experiment, there must be a coherent view of how quantities are 
defined (e.g. what precisely is a "jet"?) in order to arrive at definitive tests. 

Recent studies with high statistics of hadronic decays of the Z° from e + e~ production 



have yielded impressive new confirmations of the theory |14[ . QCD tests in pp collisions 
are not as direct as those in e + e~ owing partly to the complications associated with 
partons in the initial state and the beam fragments (the so-called "underlying event"). 
Accompanying this complexity is however a richness which allows one to attack a given 
problem in a number of complementary ways. For example, knowledge of the partonic 
density as a function of the proton momentum fraction introduces uncertainties in the 
predictive power of the theory. On the other hand the same feature allows one to obtain 
data which span a wide range of center-of-mass energies in the parton-parton frame for 
a fixed set of beam conditions. The variety and diversity of hard processes accessible in 
hadronic collisions, together with the enormous cross sections and energy reach, provide 
us with a multitude of phenomena inaccessible to current e + e~ experiments fT5fl-|20||. 

A major ingredient for the prediction of cross sections in pp collisions is the distribu- 
tion of partons inside the proton [fZT], Recently there has been significant progress in 
theory and experiment, leading to expanded measurements of parton distribution func- 
tions (PDF). Concurrently, next-to- leading order (NLO) perturbative calculations have 
been performed for most interesting processes. As a result we have collected a substantial 
body of evidence demonstrating that QCD properly describes this physics both qualitati- 
vely and quantitatively. However, there are some outstanding questions which need to 
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be resolved. The aim of this review is to present this evidence in a critical way, pointing 
out where theoretical and experimental improvements are expected or desired and where 
one can rely on QCD to extract new information. 

We review the last analyses from the CERN Collider experiments (UA1 and UA2) 
which have completed operations (1991) along with data from the Collider Detector at 
Fermilab (CDF) dating mostly from the 1988-89 running period of the Fermilab Tevatron. 
A new cycle of data taking has started in the Summer 1992 at the Fermilab Collider, with 
the presence of a new experiment, DO. Several of the results presented here are still on 
their way to press and only available in preprint form. We felt it necessary to include 
them because they often add important new contributions to the overall picture. We 
regret that the analyses from the latest set of data collected since the Summer 1992 are 
still premature to appear in this review, and we look forward to their completion for the 
important implications they will have on the study of QCD in hadronic collisions. 



2 QCD IN HADRONIC COLLISIONS 

One of the fundamental properties of QCD is the shrinking of the coupling constant as 
the energy of the interaction grows (asymptotic freedom). This implies that perturbative 
techniques can be used to study high energy hadronic phenomena. In spite of this, we 
cannot fully rely on perturbation theory (PT) because the fundamental particles whose 
interactions become weak at high energy are deeply bound inside the hadrons we use 
as beams, targets, or as observables. The solution is given by factorization theorems 



n| , 0) whereby cross sections can be expressed as the product of factors, each one 
involving phenomena appearing at different energy scales. 

In the case of hadronic collisions, the separation of the initial state evolution from the 
hard perturbative interaction can be represented, for a generic inclusive process A + B — > 
C + X, as: 

cta+b^c+x =YlJ dx 1 dx 2 f l A (x 1 , /j f ) ff(x 2 , /J, F )^ij^c(s, fx F , fx R ,a s (fx R )) , (1) 

with s = X1X2S. i and j are indices for any pair of partons (quarks or gluons) contributing 
to the process, ff~(x,fi F ) represents the number density of partons of type i carrying a 
fraction x of the momentum of the parent hadron A, and a is the cross section for the 
elementary hard process, calculable in perturbation theory []. Furthermore, the parton 
distribution functions (PDF) f(x,fi F ) are independent of the specific reaction. The uni- 
versality of PDFs is a key property, since they are not calculated from first principles as 
they contain non-perturbative information. They can then be extracted from one process 
and applied to predict rates for another one. 

The scale /zr introduced above is the scale at which the ultraviolet singularities of 
the theory are subtracted ( "renormalization prescription"), determining the "running" of 
a s (/iR,). The energy scale n F represents the freedom given by the factorization theorem to 

4 If C were a specific hadron, an independent factorization theorem for fragmentation would apply, a 
will then be the convolution of a purely partonic process with a fragmentation function describing the 
transition of a final state parton into C . 
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absorb as much or as little of the radiation from the evolution of the initial state parton 
into the PDF, including the rest in a ("factorization prescription"). 

The final result should not depend on the choice of /ip and /xr. This is the case if 
we evaluate <r, a s (/iR,) and f(x, /if) exactly. Any fixed order perturbative approximation 
will leave a residual dependence on fip and //r. This dependence is logarithmic, and 
the sensitivity of a fixed-order cross section to variations of /iF,R is usually taken as an 
estimate of the importance of neglected higher order terms. Since the two scales have 
different origins, they do not have to be the same. Nevertheless it is customary to take 
them equal and of the order of the energy scale of the hard subprocess, to avoid the 
appearance of logarithms of large ratios in the perturbative expansion and minimizing 
the effect of higher order terms. The invariance of the results under changes of /xf allows 
to formulate an equation (the Altarelli-Parisi (AP) equation, [P^]), which "evolves" the 
PDF from one scale fi-p to another. With this equation, measurements of the PDF carried 
out at relatively low values of fi in Deep Inelastic Scattering (DIS) experiments can be 
used to extrapolate the values of parton densities to the values of /i found in hard hadronic 
collisions. 

The cross sections derived at leading order (LO) have large uncertainties associated 
with the choice of /x, since the matrix elements at this order do not contain any initial 
state radiative process, and are thus independent of /ip- A dependence on /z F appears 
inside a only at the next order in PT (via the subtraction of the initial state collinear 
singularities) and a partial cancellation between /(x,/if) and ct(/zf) takes place. It is 
therefore important to have available at least the NLO matrix elements to carry out 
quantitative tests of QCD. In spite of the technical difficulties, the calculation of most of 
the interesting processes has been completed today at NLO accuracy, and new techniques 
are being developed to enable the calculation of yet higher order corrections []T3L p5| . 

Likewise, analyses of the PDFs have been carried out in recent years with similar 
precision [26]- [35] , providing the necessary elements for consistent NLO calculations. We 



refer to Ref. [pi] , |2^[ for a review of PDF measurements and parameterizations, and 
limit ourselves to point out the existence of recent data extending the measurements of 
F2,z{ x i Q) down to x — 0.008 for Q 2 as large as 5 GeV 2 [BB1 I37J. These data show clear 



discrepancies with previous extrapolations of F2 to small-x, indicating a violation of the 
light flavour symmetry in the sea densities. Nevertheless the measurements confirm [ |3"7[] 
earlier estimates of the behaviour of the gluon density, whose extrapolation to small-x 
is responsible for systematic uncertainties in the calculation of most hadronic processes. 
New fits to these data have appeared [[34], and show that the gluon density is now 
under a rather solid control in the region x > 0.01 and Q > 10 GeV (Figure [l]). This is 
the region of sensitivity for most QCD processes probed by current hadron colliders. 

The formalism described so far only allows the calculation of inclusive quantities. This 
is the case of jet distributions and correlations, or electroweak boson cross sections. The 
inclusive nature of the PDFs by itself prevents predictions on the structure of the radiation 
emitted during the initial state evolution. A more exclusive picture of the event structure 
is often required, both as a tool to understand the experimental systematics (calorimeter 
energy response, effect of particle isolation requirements, etc.), and as a way of probing 
more specific predictions of QCD (e.g. jet fragmentation properties). 

To complete the description of the event structure, a complementary approach, known 
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as shower Monte Carlo (MC), has been developed |jl0 |. Here the partons from a hard col- 
lision evolve via gluon and quark radiation, until a small virtuality scale Qq is reached 
where ct s {Qo) is large. Here confinement forces take over, hadronizing the colored partons. 
Descriptions of the hadron-formation phase can be included [^8[ |39j and tuned using a 
reference process (e.g. jet production in e + e~ collisions). These non-perturbative effects 
are believed to be universal, namely they do not depend on the hard process. The main 
features of the final state of different processes are thus accounted for by the QCD evolu- 
tion, as the distributions of the hadrons are expected to closely mimic those of the partons 
they originated from (Local Parton-Hadron Duality, LPHD 



oi 0). 

Such models have been implemented in several computer programs |4"2|-|48||. They 
differ from one another in several aspects, ranging from the accuracy of the perturbative 
evolution, to the hadronization scheme. The spectrum of the radiation is given by pertur- 
bative QCD, and, in some cases [48|, it includes all orders of leading and large classes of 



subleading soft and collinear logarithms || ]41|. MCs such as PYTHIA and HERWIG 
48|j , finally, succeed in describing typical quantum mechanical effects due to color inter- 
ference via a simple "angular ordering" prescription, which limits the kinematical phase 
space available for the emission of soft gluons from colored currents |T0, [40 . 



3 JET PRODUCTION 

3.1 Inclusive Jet Production 

The precision of QCD tests involving jets has been limited by the necessary correspon- 
dence between the final state sprays of hadrons and the partons from a hard scattering, 
whose cross sections are perturbatively calculable. There have been substantial develop- 
ments resulting from the higher center-of-mass energies and an improved understanding 
of experimental systematics. The recent calculation of jet cross sections beyond LO in 
PT has reduced the theoretical uncertainties greatly and predicts new quantities. 

The inclusive jet cross section, a(pp — * JET+X) is the most straightforward quantity 
to test. At LO (0(a*)), eight dia grams contribute to the scattering and give rise to 
two parton final states p0"| . At C(a^), the partonic cross section is directly equated to 
the measured jet cross sections. For a fixed pp center-of-mass energy, the inclusive cross 
section is a non-trivial function of two variables: r], the jet pseudorapidity (= log cot 9/2, 
where 9 is the polar angle), and the transverse energy, E t . As will be discussed later, 
the issue of how precisely one defines jet E t is important to the overall consistency of the 
comparison between theory and experiment. For now it can be taken to be the sum of the 
transverse energies of discrete sub-units, be it particles or calorimeter towers. The most 
common representation of the data is typically in terms of the differential cross section, 
da/dE t ; this is really an average of the inclusive cross section over some pseudorapidity 
interval in a detector: 

da 1 r+v da , 

dE t At]J- v dE t dr] ' v ; 

In some cases this is expressed as the cross section evaluated at rj — 0: J^l^o, assuming 
that the rate is constant in i] over a large enough interval. Most collider experiments 
report measurements in roughly the central two units of pseudorapidity. 
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As mentioned in Section 2, large uncertainties are associated in LO to changes in the 
factorization/renormalization scale \x. For a range of 2E t > /i > E t /2, the LO cross 
sections for da/dE t vary by approximately 50 %. This uncertainty is roughly a constant 
multiplier of the cross section for different Ef's , with only a modest dependence of the 
shape of the cross section as a function of jet E t . 

Despite the large uncertainties, if one chooses a renormalization scale (/i = Et/2), and 
compares data to QCD for Tevatron (y/s = 1.8 TeV)|B§, SppS (y/s = 630 GeV) [0, 
and ISR (y/s = 63 GeV) one finds an impressive agreement between the experimental 
results and the theory with only one floating parameter. This is shown in Figure p|. 

The UA2 collaboration also measured the jet cross section for different pseudorapidity 
intervals [|5l|]. Although not shown here, the agreement is good in the central region 



(| 77 1 < 0.8), but for larger values of pseudorapidity (1.2 < \r)\ < 2.0) is marginal |5T 
There is no clear explanation for this f5~lfl . 

To calculate the cross section at NLO, one must combine graphs where a parton is 
radiated and loop diagrams (Figure ^j. At this order, factors of log(/i) appear which 
cancel some of the fi dependence in a s (fi) and the PDFs. The evaluation of the full NLO 
matrix elements was initiated in Ref. |55| and later completed by Ellis and Sexton [[SlJ 



in 1986. A confirmation of these results using a different approach has come recently 



from Ref. |pq| . These works did not include an explicit calculation of the cross section. 
Whereas at LO a direct correspondence is made between jet and partonic cross sections, 
the situation is not as straightforward at NLO. In order to evaluate the cross section at 
0(ctg) (in fact even to obtain finite results f57j), one must specify what a jet is at the 



partonic level. If two partons are close together they may be merged into a single "jet" . 
Here one speaks only of jet, as opposed to partonic cross sections at both the theoretical 
as well as experimental levels. Ideally, the theoretical jet definition should thus be as close 
as possible to the experimental jet definition. 



Aversa, Chiappetta, Greco and Guillet |58| and Ellis, Kunszt and Soper |59| have 
used the matrix elements of Ref. |5(J to derive jet cross sections at 0(a^). Although 
both groups employ different computational techniques, the results have been shown to 
be numerically identical | 6"0fl . After the imposition of a jet definition (see below), the NLO 
cross sections show substantially smaller sensitivities to renormalization scale variations 
than at LO. Over a range of renormalization scales close to the hard scattering scale 
(E t /4 < fi < Et), the uncertainties in the cross section have been reduced from 50 % 
to 10 % over most of the range of accessible E t . The inference is that the effects of still 
higher order contributions are rather small at 0(a z s ). Figure f| shows the variation of 
the 0(0*1) and 0(a 2 s ) cross sections with \x for 100 GeV jets. At LO, one finds a large 
monotonic variation of the cross section with fi, whereas at 0(a1), the negative log(/i) 
contributions from the virtual terms reduce the cross section at very small fi. Note that 
the sensitivity of the cross section at renormalization scales near the hard scattering scale 
is greatly reduced at NLO. 

Several experimental jet algorithms have been employed. When cross sections are 
derived only at LO and when uncertainties are large, these differences can be forgiven; 
even so, comparisons between experiments are rendered difficult. For example, UAl 
and CDF employed cone algorithms, whereas UA2 used initially a nearest neighbor 
algorithm f6~2fl . A typical hadron collider algorithm is the "cone" algorithm, which has 
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been suggested as a standard for pp experiments [p^] . It operates in a space defined 
by pseudorapidity and azimuth (rj — 0) on particles, or partons or calorimeter towers, 
depending on the specific application. In this metric, one can define a jet to be the 
partons or particles found in cones or, more precisely, circles of radius AR = \J A</> 2 + Ar/ 2 . 
The transverse energy, E t , is the sum of the transverse energies of particles, partons or 
calorimeter towers inside a fixed radius. The direction of a jet in rj and can be defined 
as the E t weighted centroids: 

E>« = (3) 

i 

<t> jet = CL E M/El et (4) 

i 

Viet = C£EuV l )/Er t (5) 

i 

The above description is not complete, however. It does not tell where to initially 
place the cones to form the above quantities, and does not describe how to handle cases 
where cones overlap ("merging"). In the case of experiments employing calorimeters, 
the initial jet direction can be defined by towers with E t above a given threshold (seed 
towers). An iterative approach can be adopted to find a stable center of the cluster by 
successively recomputing the cluster centroid until the list of towers or particles in the 
cone is stable [p3|| . If two jets are greater than one cone radius apart, but less than two 
radii {i.e. R Q < Ai?i j2 < 2R Q ) should they be identified as one or two jets? The inherent 
difference between two partons in a calculation and calorimeter towers in an experiment 
can make it difficult to achieve a precise uniformity in the jet definition. 

Both CDF and the UA2 experiments have measured jet cross sections using cone algo- 
rithms with R = 0.7 and 1.3 respectively, with reduced uncertainties p5| , [fjl] . Although 



the UA2 results are not shown, there is good agreement with the 0(a^) predictions de- 
spite the fact that the calculations do not apply strictly for R > 7r/3 [|59|] . The dominant 
experimental uncertainties are associated with the hadronic energy scale. The calorimeter 
response to jets, particularly the hadronic component, is difficult to calibrate in an abso- 
lute way. There are no test beams with monoenergetic sources of jets, so the calorimeter 
response must be derived from a convolution of the calorimeter response to hadrons of 
varying energy (including 7r°'s) with the jet fragmentation spectrum. Although the re- 
sponse can be checked with sources such as jets recoiling against direct photons, there 
is no unimpeachable source on which to calibrate. Typical energy scale uncertainties are 
now ^4-15% in 5E/E. A systematic shift in energy scale is equivalent to an uncertainty 
in the cross section. Since the cross section is typically a steeply falling function of E t , 
following a power law spectrum of E^ 5 , the resulting uncertainty in jet cross section is 



20-75%. Recent work by both the UA2 [51] and CDF [65] collaborations pressed the 



lower bounds of these uncertainties, thus improving the level of comparison to theory. 
The uncertainty can be expressed as an overall multiplicative factor which is independent 
of jet E t (20% and 35% for CDF and UA2 respectively) and a smaller term which is E t 



dependent and can be roughly 5% |65| 



The agreement appears to be very good on a logarithmic scale. To illustrate significant 
features of the comparison, however, one can plot the cross section on a linear scale, as a 
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ratio of (Data - Theory) /Theory as a function of jet E t . Figure [| shows such a comparison 
for CDF data f6|. The Q CD °( a s) prediction for // = E t /2 is defined to be the "Theory" 
or on this plot for the purposes of normalization. The data have uncertainties factored 
into a combination of the E t dependent systematic and statistical uncertainties which are 
displayed on the error bars, and an E t independent component which is 20 %. Figure [5] 
(a) illustrates the improvement in the uncertainty associated with theory for a variation 
of E t /A < fi < E t . One can see that the uncertainties are substantially reduced at 0(af). 
Figure |] (b) shows the effect of different PDFs |31], [£J on the predicted cross section. As 
one can see there is some dependence on the shape of the derived cross section on the 
choice of PDF, however, the overall agreement is quite good. 

This does not exhaust comparisons at 0(al). Figure [| shows the variation of the cross 
section with cone size for E t =100 GeV jets from CDF compared with 0(a^) predictions 



[66], |68| . The data display the statistical errors only, but the ps 23 % systematic 
uncertainties are largely independent of R. Since there are only two back-to-back partons 
in the LO calculation, one can only predict such a variation beginning at NLO. An 
interesting feature of the calculation is the minimal sensitivity to [i for a cone size of 
R = 0.7, whereas the sensitivity is much greater both for R < 0.5 and R > 0.9. From this 
standpoint R = 0.7 represents an "optimal" cone size for comparison to 0(a^) predictions. 
The data appear to be in rough agreement with at least one of the QCD predictions 
(/x = E t /4), but on the whole, there seems to be a trend for the data to show a slightly 
steeper dependence on R than the theory predicts. 

A quantity related to the variation of cross section with cone size is the jet profile. 
To measure this, one can pick a large radius (R = 1.0), and then examine the fraction 
of the jet E t contained in a smaller sub-cone of radius r: E(r, R, E t ). CDF measured 
this quantity using charged particle tracking data because it is more fine-grained than 
calorimetric information. Figure |7| shows a plot of F(r, R, E t ) from CDF data. Also 
shown are the predictions of 0(a^) QCD [67] for different choices of renormalization 



scale. It is perhaps surprising that the data are so well described at the level of just 
one gluon bremsstrahlung when there are typically 10 charged hadrons in a typical jet. 
Since 0(a^) is the lowest order at which one can speak of a jet profile, the sensitivity to 
renormalization scale is fairly large. 

There is an apparent contradiction between the profile measurement and the variation 
of the cross section with cone size. One naively might expect that if there were good 
agreement between data and theory for one quantity, having chosen a renormalization 
scale, that there would be a good agreement for the other. This expectation is based 
on the assumption that the variation of cross section with cone size just depends on the 
energy flow within the cone. This assumption is not valid, however. In the inclusive 
measurement, jets are clustered independently for each cone size, R, chosen, whereas 
for the jet profile, only a single cone of R = 1.0 is used. The main difference is the 
"merging" step. Ellis, Kunszt and Soper |67] have examined the effect of merging in the 
0(al) predictions. As discussed above, there is an ill defined region where two partons 
may be separated by a distance R a < AR < 2R Q . In order to mimic the merging in the 
experimental algorithm, partons are merged into a single jet if they have a separation 
AR < R sep . The calculation implicitly had R sep = 2.0. However as seen in Figures ^| 
and [7| a value of R sep = 1.3 and choice of fi = E t /4 fit both distributions |)7]]. Although 
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one has obtained consistent results, one has done so at the expense of adding a tunable 
parameter to the theory. 

The ambiguities related to the prescription for the merging of jets are absent in the 
class of jet definitions generally used by the e + e~ experiments. The prototype of these 
jet definitions is provided by the JADE algorithm p5[, which builds clusters of charged 



particles according to an invariant mass cut. The invariant mass normalized by the 
center-of-mass energy, = Mfj/E^ m , is used to define jets as distinct objects. My 
is the invariant mass of pairs of particles or of a particle and a cluster. At each step 
of an iterative procedure, the pair with the smallest yij is merged into a new cluster if 
Vij < Vent- If no pair is left passing the cut, all remaining clusters are called jets. The 
leading weakness of the JADE algorithm from the point of view of pp, pp and ep colliders 
is that all particles are associated to some jet, including those coming from the underlying 
event and which do not belong to the hard process. 

Improved versions of the JADE algorithm have recently been proposed |T0fl , which 
reduce the sensitivity to the jet definition under hadronization corrections, and make 
it possible to resum large classes of leading and subleading perturbative corrections in 



the theoretical calculations. These prescriptions can be extended |7T| to processes with 
hadronic the initial states. In this formulation they provide an unambiguous prescription 
for the merging of jets and allow the universal factorization of initial state collinear sin- 
gularities, minimizing the contamination from the hadron remnants and the underlying 
event. The similarity with the e + e~ jet definitions will make it possible to compare jet 
properties between e + e~ and hadron colliders in a consistent and universal fashion. No 
complete phenomenological study of this new algorithm is available as yet, but we hope 
that progress will be made soon (S Ellis, Z Kunszt & D Soper, personal communication) 
and that experimental measurements will follow as well. 

3.2 Xt '. Jet Scaling with s 

If one plots the inclusive jet cross section in terms of two dimensionless variables, the 
"scaling" hypothesis predicts an independence of pp center-of-mass energy, s. In reality, 
the evolution of PDFs and ot s with the hard scattering energy scale causes a violation of 
scaling for the inclusive jet cross section. 

To test scaling, one typically plots Ef times the invariant cross section (Eda 3 /dp 3 ) 
as function of x t = 2E t / y/s to obtain two dimensionless quantities to express the jet 
cross section. If scaling were valid, cross sections measured in this way at any y/s would 
all fall on a single universal curve. QCD, on the other hand, lifts the degeneracy. The 
predicted ratio of cross sections at two different center-of-mass energies as a function of 
Xt is relatively insensitive to choice of PDF, renormalization scale or the order of the 
calculation, making it a relatively solid test of the theory. Independent measurements 
made at the SppS, ISR and Tevatron showed rough agreement with QCD scale breaking 
5T|j . CDF |72[ have recently compared jet cross sections at y/s = 546 and 1800 GeV as 



a test of x t - When the cross sections are measured in one experiment, a large part of 
systematic uncertainties (e.g. hadronic energy scale) cancel when the ratio of the cross 
sections is taken, improving substantially the level of comparison. 

Figure [8] shows the ratio of scaled cross sections as a function of x t for CDF data taken 
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at both center-of-mass energies. The error bars show statistical uncertainties, and the 
shaded area indicates an overall systematic uncertainty in the ratio. The data are clearly 
inconsistent with scaling (Ratio=l). The data do not exhibit a wonderful agreement with 
QCD either. 

The discrepancy with QCD is about 2 standard deviations in the systematic uncer- 
tainty, which is not sufficient to indict QCD by any means, but is curious. At the moment, 
there is no obvious explanation for such a discrepancy. PDFs in the relevant x range 
(x > 0.1) have been measured by a number of DIS. Further running of the Tevatron 
collider at a lower center-of-mass energy could shed light on this. 



3.3 Two-jet Distributions 

The 0(a^) predictions have been extended to measurements where one defines a two-jet 
inclusive final state. Since it is impossible to either measure or calculate states beyond LO 
with two and only two jets, due to soft radiation, one can form quantities from the leading 
two jets, and ignore other energy in the event. Recently, Ellis, Kunszt and Soper |6(| have 



extended the 0(a^) calculations to predict the two jet invariant mass and center-of-mass 
angular distributions. Such distributions are sensitive to the presence of deviations from 
QCD arising from quark compositeness, technicolor [73f and axigluons [[74[ . 



The CDF two-jet invariant mass distribution Mjj is defined as: 



Mjj = v (-^i + E 2 ) 2 — (pi + P2) 2 (6) 

where E{ and pi are the energies and momenta of jets 1 and 2. Note that the effective 
masses of the jets enter into the determination of Mjj. The jet mass, an internal quantity, 
can be associated with gluon bremsstrahlung within the clustering cone. The CDF Mjj 



cross sections were determined for cone sizes of 0.7 and 1.0 ||7q| . For the HMRSB [[IT 
and MT SI [33] PDFs the 0(a^) predictions appear to be in good agreement with the 



data for a clustering cone of 1.0, but disagree at about the 2 standard deviation level for 
a cone of 0.7 [75]. 



The dijet angular distribution has likewise been calculated at 0(oP s ) Since in- 

variant mass and cos(6 1 *) are independent variables, the data can be placed in different 
bins of Mjj. Here 6* is the center-of-mass polar scattering angle. Since the cross section 
is dominated by t channel exchange, it rises very rapidly with increasing cos(6*) and it is 
more convenient to plot the data as a function of the variable Xi defined as: 

= i + cose* 

X ~l-cos9* [) 

If plotted versus Xi the Rutherford scattering pole is taken out. There is a rise in cross 
section for x ~ 1 (90 degrees) associated with the contribution of s-channel scattering. 
Figure 19] shows the results of an analysis of dN/dx by CDF using a cone size of 0.7 [75 



The data are divided into three bins of Mjj. One can see that the data are well described 
by both 0(a 2 s ) and 0(af) QCD predictions. The data are separately normalized for each 
bin of Mjj. 

The effects of quark compositeness would be to increase the amount of data found 
near x ~ 1- Quark compositeness is typically parameterized in terms of a four-fermion 
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interaction with a coupling inversely proportional to a characteristic energy scale (related 
to the "size" of the quark), A c [[77], ^5[. Such an interaction gives rise to an isotropic 
distribution in the center-of-mass system, and also contributes a rising cross section at 
large E t or Mjj. In order to search for compositeness, one can take several bins of Mjj 
and examine the dijet angular distribution in each. Compositeness could be manifest as 
an increase in 90 degree scattering in the highest Mjj region, while the remaining data 
should be well described by QCD. The CDF data have allowed limits to be placed on 
A c > 1.4 TeV using the inclusive jet data p5| , and at A c > 1.0 TeV using the angular 
distribution If76 



3.4 Jet Fragmentation 

As discussed in Section 2, some aspects of jet fragmentation reveal the underlying QCD 
mechanisms, particularly when one assumes that the behavior of hadrons in jets mimics 
features of the partonic emission. On the basis of this, one expects for example that 
gluon jets will have softer fragmentation than quark jets and that average multiplicities 
will increase with energy. 

Studies of jet fragmentation in hadronic collisions have been performed in the past 
by UAl and UA2[j7lJ, providing the first indications that jets in pp reactions have higher 
multiplicities than in e + e~ annihilation. 

The most notable quantity to study is the jet fragmentation function, which describes 
the probability of finding a hadron carrying a given fraction of the jet's momentum. This 
is typically described in terms of the component of hadronic momentum parallel to the jet 
axis: p\\. The charged particle fragmentation function, F(z) , is defined as (z =pn/pj et ): 

„, . 1 dN c h , s 

The evolution of F(z) as a function of the hardness of the primary collision is a good 
test of QCD, and, in principle, can be used to extract a s . The same mechanism for the 
evolution of parton densities, namely soft and collinear parton emission, is responsible for 
the logarithmic evolution of F(z). Figure |TT| shows the evolution of different bins of F(z) 
as a function of dijet invariant mass {Mjj) from CDF data JFjJ. The data agree well with 



a logarithmic evolution with Mjj and have a distinct similarity with data from e + e~ |80 
which are plotted as a function of y/s. Mjj appears to be a sensible variable to express 
this evolution insomuch as it is a measure of the hardness of the scattering, particularly 
in the central pseudorapidity region. 

Notice however that the e + e~ and pp curves do not match: the e + e~ curve correspond- 
ing to the lowest z bin extrapolates below the equivalent CDF curve. This behaviour is 
consistent with the notion that jets in hadronic collisions are mostly produced by gluons, 
while in e + e~ they come from the evolution of quarks. We should however point out that, 
as noticed in a previous section, jets are defined according to different algorithms in e + e~ 
and hadronic collisions, and unless a common definition is provided it is not possible 
to draw quantitative conclusions from these comparisons. Nevertheless it is encourag- 
ing that, as shown in Ref. J79|, calculations based on the HERWIG Monte-Carlo are in 



agreement with the measured inclusive fragmentation function. 
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3.5 Multijet Final States 

In all of the above, one has considered final states where two jets predominate. Predictions 
for multijet final states are obtained either from QCD shower MC programs or from fixed 
order PT. For the latter, one can obtain finite cross sections by limiting the minimum E t 
and opening angle of partons in order to stay away from soft and collinear singularities. 
For states selected where the partons are stiff and widely separated, one expects that tree 
level predictions should be reasonably faithful. There is no NLO calculation available 
for Nj et > 2, therefore since one does not have the log(/i) cancellation that appears at 
NLO, there is a substantial sensitivity in the predicted cross sections to a variation of the 
renormalization scale. This is because the cross section is of order a s (fi) N where N is the 
number of final state partons. Any uncertainty in the scale will hence be multiplied by a 
large factor in deriving cross sections. 

The tree level matrix elements commonly in use are based on calculations by several 
groups^, and have been included in numerical programs [[8^, 83] which are currently used 



by the experiments. Because of the complexity of the results, techniques have been 



developed to provide reliable approximations to these matrix elements ||34|-|[88||. The 
testing of these approximations using current data is very important, because rates for 
multi-jet production at the future colliders will be extremely large and fast but reliable 
numerical simulations will be required to evaluate them. 

To start with, the topologies of multijet final states appear to be well predicted by 
the tree level calculations. Several examples can be seen in both CDF and UA2 data. 
CDF examined the topology of three jet events with high statistics and in regions of 
uniform acceptance. They found a very good agreement with tree level predictions |63 



In particular, there is a distinct difference expected between three jet topologies for events 
initiated by gluon-gluon and gluon-quark collisions versus those from quark-antiquark 
collisions. The data appear to be in good agreement with the expectation that most of 
the three jet final states come from gluon-gluon and gluon-quark collisions [p3 |. 



Both UA2 and CDF explored in some detail the structure of four jet final states 



89|. p0|| . These studies are partly motivated by a search for double-parton processes where 



two uncorrelated 2 — > 2 scatters occur, producing four jets in the final state. Figure [12 



shows the angular separation of all pairings of jets from four jet events in CDF data 



90[| and compared with the results of predictions based on the exact tree-level matrix 



elements p2| |, which reproduce the data very well. In contrast to QCD production where 
the four jets have no intrinsic correlation, the double-parton events are expected to have 
jet pairs which approximately balance in transverse momentum. One expects that the 
cross section for double parton events would scale like the square of the dijet (2 — > 2) 
cross section, normalized by a factor that is comparable to the inelastic cross section 



a DP = adijet X adijet (9) 

2(y eff 



where a e ff is expected to be roughly 10 mb ||90|| . The factor of 2 in the denominator 
is included to account for the Poisson nature of chance of a double parton interaction 
The Axial Field Spectrometer collaboration reported a significant double-parton 



'For a review of these techniques and for a complete set of references, see Ref. fil 
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cross section, with er e // = 5 mb [I5I|, whereas the UA2 collaboration did not find any 
evidence for the process and set a limit of a e ff > 8.3 mb (95 % C.L.) |j89|pt Note that 
since a e ff appears in the denominator of equation (Q), a larger number implies a smaller 
(Jdp- Finally, CDF reported an effect at roughly a 2.5 standard deviation level, with 
a e jf = 12.ltl°4 mb. It is estimated for four jet final states at the SSC, that double 



parton scattering will dominate for jet p^s less than 40 GeV ||90|| . If, in the future, a 



sizable effect is observed, it is possible one may obtain unique information on correlations 
between partons in the proton from double parton scattering. 



The UA2 collaboration studied the cross section of events with up to six jets. Figure 13 



shows the jet p t distributions for 4, 5 and 6 jet final states, compared with various tree level 
predictions. Notice that the normalization of the theory curves is absolute. Considering 
the complexity involved in these calculations, the agreement with data is remarkable and 
extremely encouraging in view of the potential applications of these calculations to the 
study of multi-jet phenomena at the future hadronic colliders LHC and SSC. 

As an alternative to exact tree level calculations, and to get a more exclusive descrip- 
tion of the events, one may employ shower MC's. In this approach multiple jets can 
appear when branchings with large transverse momentum relative to the leading partons 
take place. Given the approximations involved in evaluating these large p t branchings 
inside the MC, an assessment of the reliability of this approach is in order. 

CDF recently performed a detailed comparison of the characteristics of events with 
high total transverse energy with the HERWIG event generator combined with a realistic 



detector simulation ||93|| . The events were selected by requiring that the total transverse 
energy be in excess of 400 GeV. Events with up to 6 jets were observed, and the HERWIG 
generator does an impressive job in reproducing a very large number of distributions, 
such as the jet multiplicity as a function of different jet-p t thresholds (Figure the 
jet profiles, invariant masses of various combinations of jets. Such studies illustrate the 
power and accuracy of event generators to reproduce event characteristics. 

Although they agree in rough detail, there are some significant differences among 
some of the MC event generators. One of the most relevant differences is how color flow 
is handled. In particular, the dynamics of color flow leads to the need of the angular 



ordering of QCD radiation in shower MC's |10|- The emission of radiation is related 
to the color structure of the hard scattering process. Some event generators, such as 
ISAJET ignore the connection between radiation and the hard scattering, whereas other 
generators, such as HERWIG have explicitly built in the color flow connection to better 
reproduce event properties. The differences between the coherent and incoherent emission 



has been studied extensively in e + e JT4|], but not as much in pp collisions. This is largely 



due to the inherent problems in distinguishing the soft particle flow associated with the 
hard scattering from that associated with the underlying event. For sufficiently large 
momentum transfers, however, one expects the radiation effects to become visible as jets, 
which are more readily associated with the hard scattering process. CDF studied the 
angular distribution of the third highest E t jets in events with two high E t leading jets 



95]. These studies show significant differences between the predictions for ISAJET and 



6 The findings by UA2 and AFS are not necessarily inconsistent: on one side the x-range probed by 
the multi-jet configurations at the two energies of 63 and 630 GeV is very different. On the other, at the 
time of the AFS analysis the exact predictions for the QCD four-jet production were not available. 
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HERWIG, where the data are in much better agreement with HERWIG, indicating that 
for some measurements, color coherence effects cannot be neglected. 



4 HEAVY FLAVOR PRODUCTION 



Heavy quark production in high energy hadronic collisions constitutes a fundamental arena 
for the study of perturbative QCD. Of particular importance is the role played by mq. 
Only in b quark production does one have today the unique situation where rrig 3> Aqcd- 
The prediction of heavy quark production cross sections in hadronic collisions has far 
reaching implications. Discovery reaches and limits for the "top" quark depend on reliable 
estimates from PT. The observability of CP violation in B mesons at hadron colliders 
depends, to a large extent, on the production cross section and correlations between 
the B and B. Recent years have witnessed remarkable progress both in the theoretical 
understanding of the production mechanisms [0 and in the experimental capability to 
probe them via independent and complementary observations |98] . 

The mass of the heavy quark Q provides a natural infrared cutoff in the evaluation of 
the production rates and multiplicities. Complete NLO calculations are available today 
for the total [pj , one-particle-inclusive ||100|| and two-particle- inclusive ||101|| cross sections. 
Production of heavy quarks in the perturbative evolution of high energy jets has also been 
studied, and LO expressions for the heavy quark multiplicities are known ||102| . 

The non-perturbative corrections which are required to derive the production proper- 
ties of observable heavy flavored hadrons hq are suppressed by powers of Aqcd/wq- For 
production at large p t , the factorization theorem guarantees the existence of a fragmen- 
tation function D® Q (z,fi) which models the fraction of momentum of the heavy quark 
retained by the heavy hadron: 



d 3 p 



E'd 3 



d 3 p' 



—D? 



, . dz 



(10) 



where p = zp' and ctq is the elementary cross section for the production of the heavy quark 
Q, calculable as a perturbative expansion in a 8 . The evolution of the fragmentation func- 
tion with the factorization scale \l obeys the AP equation 
which is given by D?(z,rriQ) = 5(1 — z), up to non-perturbative effects 



with a boundary condition 



These 

non-perturbative effects obey a scaling law in itlq and can therefore be parametrized in 
a phenomenological way by fitting, e.g., e + e~ data ||104| , |105|| . With this additional in- 
put, non-perturbative corrections to Equation [l^ are suppressed by powers of p t . The 
evolution of D® Q (z, /i) with \i is known today up to NLO in PT [ |106| . 

When applied to the energy of the current hadron colliders, these results are believed 
to provide a reliable description of the production properties of very massive quarks - 
e.g. the yet undetected top. In the case of charm and bottom, the situation is more 
delicate. In fact production of c and b quarks is dominated by gluon fusion processes 
(99 ~^ QQ) an d the distribution of gluons inside the proton is probed at values of x close 
to the boundary of current DIS measurements. Furthermore the NLO contribution is 
larger than the LO result, and very sensitive to the input scale fi. Significant corrections 
are thus expected from yet higher order terms. These corrections arise from a class of 
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diagrams with t-channel gluon exchange first appearing at NLO ||107|| (Figure [T3]) . They 
lead to terms proportional to powers of a s log(s/mg 2 ) | 99| , which might dominate at higher 



energies, as well as becoming non negligible in the case of top production at supercollider 
energies. Techniques exist to resum these large logarithms ||108|| , and have been extended 



for application to this specific problem [ 1TDDL |110| , |111[ | . Comparing the results of the NLO 
predictions with the available data and verifying whether the resummed calculations can 
explain possible differences is therefore of utmost importance as a test of QCD per se and 
as a milestone before extrapolation to higher energies. 

4.1 Bottom Production 

There are several different channels which allow the detection of b quarks. Fully re- 
constructed exclusive decays of 6-hadrons allow the unambiguous tagging of a 6-quark, 
together with a precise measurement of the hadron momentum. Viable examples are pro- 
vided by B ± -> J/ipK* |]]| PH, B° -» J/ipK* [TrjJ and A b -> J/ipA Due to the 



small branching ratios (BR) and detection efficiencies, these channels are only accessible 
near threshold (p t = 0(mi)), where the production rate of b quarks is more abundant. 
The region of small p t is expected to be more sensitive to the uncertainties in the calcula- 
tions mentioned previously and is therefore potentially more interesting for critical tests 
of QCD. 

At larger values of pt (typically above 10-j-15 GeV) semileptonic decays become the 
leading tool to study b production. Neglecting detector backgrounds, and neglecting W, Z 
and c decays, b quarks are the most abundant source of high p t leptons. Several techniques 
can be employed to subtract the above backgrounds [p8 [. Backgrounds from Z's, W's and 



continuum Drell-Yan events can be identified because single leptons from these processes 
are more isolated than leptons from heavy quark decays, surrounded by the fragments of 
a jet. In addition, lepton pairs from Z's can be eliminated with a cut on the invariant 
mass of the lepton pair, and W's can be identified by the large transverse mass of the iv 
pair. 

For pt values larger than 10-=-15 GeV, the c and b cross sections are comparable. Since 
b quarks undergo a harder fragmentation into hadrons compared to c quarks, and since B 
hadrons have a larger phase space available for the decay, we expect the c contamination 
to contribute only a fraction of the total lepton yield. This fraction can be precisely 
estimated by studying the transverse momentum of the lepton relative to the direction of 



the jet in which it is imbedded [[112 



Furthermore, the b component can be determined by tagging charmed hadrons (say 
.D's) inside the jet and with the correct charge correlation with the lepton itself, e.g. 
e~D° as opposed to e~D°. UAl has also pioneered a technique based on the detection 



of a second lepton in the event [ 112[| . This second lepton comes either from the charm 
emitted during the decay of the b into the leading lepton, or from the decay of the second 
b in the event. In the first case we have a low mass dilepton pair, and the measured rate 
can be directly related to the b cross section. In the second case we have a high mass 
dilepton pair, and the extraction of the inclusive b cross section requires an understanding 



of the correlations between the two heavy quarks in the event |120| . 



The advent of new technologies, such as secondary vertex detectors capable of isolating 
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the charged particles coming from the displaced vertex of a B decay, will provide further 
tools to strengthen the capability of hadron collider experiments to tag b hadrons and 
study their properties. 

Unlike other inclusive measurements (e.g. direct photon, jet production), inclusive b 
cross sections are reported as a function of the integral cross section above some pt- This is 
done to minimize systematics associated with the b fragmentation and decay. The effects 
of these two effects must be unfolded in order to obtain a b cross section from the observed 
lepton spectrum. The results of the measurements by UAl [ |112|| and CDF [ |114] |-[ |TT6|j are 



collected in Figure 16]. The two solid lines represent the NLO QCD prediction ||10C 



obtained using PDF's from the most recent MRS fit (set DO) and two different values 
for /i and Aqc^I]. This band is supposed to represent an acceptable range of variation for 
the input parameters of the NLO calculation. The value of was fixed to 4.75 GeV. 
A variation of the mass in the range 4.5 GeV < < 5 GeV only affects the result by 
no more than 20% in the region p t < lOGeV^, and of the order of few % above 20 GeV. 
Two features are to be noticed. First of all, the theoretical uncertainty is rather large, 
significantly larger than the uncertainties encountered in the case of the NLO inclusive 
jet cross section. Secondly, while the UAl data fall well inside the theoretical band, the 
CDF points are systematically higher, with deviations of up to a factor of 3 for the low-p t 
points. 

No satisfactory explanation for this discrepancy is available as of today, even though 
at least two suggestions have been put forward. First of all the gluon momentum fractions 
x probed by the CDF measurements are significantly smaller than those probed by UAl. 
Attempts have been made [ |118| ] to explicitly include the CDF b data in global fits of 



the gluon density. These attempts have not led to a complete solution of the problem. 
An explanation of this can be found in the following observation | |119|| : the region in 



x which is unexplored even by the most up-to-date DIS data is x < 0.01; using the 
available extrapolations of the gluon densities below this value, the contribution to the 
cross section for 6's with p t > 10 GeV coming from the region x < 0.01 is only of the order 
of 20% (Figure [17]). Therefore only large differences in the extrapolation could explain the 
observed discrepancy, and such differences are difficult to achieve because of the global 
constraints posed by the measurements of gluon distributions at larger values of x, such 
as momentum sum rules. 

An alternative explanation could be provided by the presence of the large log(s/m 2 ) 
corrections mentioned previously. The studies in Ref's [ l(J9]-[ 111] have led to a general 



reformulation of the factorization theorem for application to processes where initial state 
gluons with small momentum fraction x are involved. The result can be expressed in terms 
of gluon distributions depending not just on x and fi, but on the transverse momentum 



k as well fi09fl : 



o (s 



rl rl roc roc 

/ dxi / dx 2 \ dk\ / rfk2^ r (xi,ki, / u)J'(x2,k2,/i)/(s, k 1; k 2 ), (11) 
Jo Jo Jo Jo 



7 It is worth pointing out that the values of Aqc d extracted from fits to DIS data are systematically 



lower than those obtained from precision measurements of jets performed at LEP [117|. The differences 
are of the order of two standard deviations. Using for a s the values extracted by LEP experiments would 
increase the predicted b cross sections by an additional 20%. 
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where the functions T describe the transverse momentum distribution of gluons with 
longitudinal momentum faction x. I, referred to in the literature as the impact factor, 
represents the gauge invariant elementary cross section for the process gg — > QQ with 
initial off-shell gluons of virtuality — k . An intuitive physical interpretation of this result 
is the following: at small x and for \x 3> ^.qcd, gluons are more likely found in a peripheral 
branch of the initial-state evolution tree. In other words, the multiplicity is dominated 
by processes where the degradation of the gluon momentum down to a fraction x took 



place via a large number of successive splittings (see Figure [18]). Since at each splitting 
the gluon acquires some transverse momentum k, k will build up during the evolution to 
small-x; for x small enough the transverse momentum will not be negligible with respect 
to the scale of the hard process, \i. Therefore the description of the gluon density at small- 
x should depend on k as well as on x and /x, and its evolution equation cannot neglect 
the transverse degree of freedom. An evolution equation for the density T[x, k, /j) can be 
formulated ||108| , extending the standard AP equation. This evolution equation resums 
the leading (a s log(s/m 2 )) n terms which appear in the perturbative expansion for the hard 
scattering cross section and allows them to be absorbed into T{x, k, /j), provided one uses 
the impact factor I rather than the standard on-shell matrix element in the expression for 



the cross sections, Equation [11]. The result of this approach cannot be simply estimated 
by varying the renormalization scale \x within some range, because the impact factor and 
the k-dependent density contain information beyond what available in the standard NLO 
calculation; this could explain why even the change of /i in the rather extreme range of 
mj/4 < fi < tut cannot reconcile the NLO prediction with the data. 

The main physical consequence of this picture is that small- x gluons involved in a 
hard scattering at a scale \x will have an intrinsic transverse momentum of the order of 
\i itself. This additional transverse momentum will smear the pt distributions obtained 
from a pure NLO calculation, but complete calculations of this effect are not yet available. 
Explicit estimates exist [ 109 ] of the corrections to the total cross section resulting from 
Equation [11]. At Tevatron energies these corrections amount to approximately 50% of 
the NLO total cross section. While this effect seems insufficient to explain the observed 
discrepancy, one should keep in mind that the smearing induced by the effective intrinsic 
p t introduced by Equation |TT] could very well push most of this contribution to values of 
Pt > nib, where the NLO cross section is only a fraction of the total. 

While we await for more calculations, it is worth exploring additional consequences 
of this scenario. In addition to pushing the measurement of tVs to smaller values of p t , 
it is useful to study correlations between the b pair. NLO calculations exist for these 
correlations |[101|| . If the small- a; effects behave as indicated, we would expect to observe 
a flattening of the A(f> and pt bb distributions relative to the NLO prediction. Here A<f) 
represents the difference in azimuth between the b and the b, and p t bb represents the 
transverse momentum of the pair. The flattening would be caused by the additional 
intrinsic p t due to the gluon transverse momentum k. 



The A0 correlations have been studied by UA1 ||120 , indicating a good agreement 
with the NLO calculation. This result does not resolve the issue, however, because the 
agreement of the NLO b cross section with the data suggests that the energy at UAl is 
below the threshold for the possible onset of these new small- x phenomena. 
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4.2 Charmonium Production 



In this Section we review the status of the measurements of production cross sections for 
charmonium resonances such as the J/ip. The theory of quarkonium production |L2j| is 
on a less solid ground than the theory of open heavy-quark production. Production cross 
sections are evaluated by convoluting the cc matrix elements with the non-relativistic 
charmonium wave function, parametrized in terms of the decay widths of the relevant 
(J, L) state. The QCD radiative corrections to the LO processes have not been evaluated 
yet. 

The observation of J/ip's is however an important ingredient in the study of b pro- 
duction. On one hand, a significant fraction of the detected J/ip's comes directly from 
6-hadron decays rather than from prompt charmonium formation \\122[ |123|1 . In fact the 
J /ip form factor inhibits production with p t 3> m c . On the other hand, 6-decay final states 
involving a J/i[) provide unique tags in the search of yet unobserved or rare 6-hadrons 
(such as B s , B c , A&) as well as in the detection of CP asymmetries (e.g. from Bd — > J/ipK® 
decays [(5(J.) A coherent picture of the production of both b and J/ip in hadronic collisions 
will therefore provide not only a significant test of QCD, but also the starting point for 
important studies of the Standard Model. 

Figure |19| shows the inclusive p t differential distribution for J/i/j's measured by UAl 
123| and CDF [124 . We superimpose the result of a QCD calculation | |119| | based on 



the LO matrix elements given in Ref. ||121|1 for the direct charmonium production, plus 
the contribution from B decays evaluated using NLO matrix elements jlOOfl , convoluted 
with a Peterson fragmentation function and the experimentally observed B — > J/ip decay 
spectrum. The theoretical error band is evaluated using the same range of parameters 
Aqcd and \i employed before in the study of the b cross sections. Notice that changing 
\i for the direct charmonium contribution causes a variation ranging from a factor of 7 to 
10, depending on p t . This indicates that the LO prediction for direct charmonium is very 
poor, and very large NLO corrections should be expected. 

In the case of UAl the data fall all inside the theoretical band, while again CDF 
shows a production rate larger than expected. A similar feature is observed in the CDF 
measurement of the ip{2S) p t distribution ||124|| . 

An important parameter is the fraction of J/^'s coming from b decays, Jb- This 
number allows to extract a b cross section from the observed J/ip production rate. Notice 
from the theoretical curves in Figure [L9| that f B is very sensitive to the parameters used 
for the evaluation of the two contributions. 

fs can be extracted experimentally, for example by separating the direct J/ip's from 
those due to B decays via the observation of the displaced vertex from which the J/ip 
originates, due to the long B lifetime. UAl measured fs (32% for p t (ip) > 5 GeV | |123| ) by 
assuming that direct J/ip's are isolated while J/^'s from B decays are not. This number 
is consistent with the estimates provided in [|119|| . 

The assumption used by UAl to extract fs might not be correct if other produc- 
tion mechanisms were responsible for direct quarkonium production, such as for example 
gluon — > J/ip fragmentation [ |127|| . It is reasonable to expect that at some value of pt the 
dominant production mechanism for charmonium states will indeed be via gluon fragmen- 
tation. The main reason being that direct production as described by the LO mechanisms 
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inhibits production at large pt via a form factor suppression (the probability that a char- 
monium bound state will hold together when produced directly in an interaction with a 
large virtuality scale is highly suppressed). The fragmentation functions for the creation 
of S-wave charmonium (rj c and J/ip) in a gluon shower have recently been calculated 
127| | and those for the creation of P-wave states (x) will soon be available (E Braaten & 



TC Yuan, personal communication). 

These calculations can be used to extract the fragmentation contribution to char- 
monium production in the regions of p t explored experimentally, and to verify whether 
this process can account for the large observed rates. The experimental detection of non- 
isolated J/^'s from a primary vertex, therefore not coming from B decays, would indicate 
that these processes are indeed present. 

Measurements of the decay-vertex position of the ip(2S) would provide evidence in 
favour or against the current belief that most of them come from B decays. If the gluon 
fragmentation mechanism were important, it would appear with a signal of non-isolated 
prompt ip{2S). 

Similarly interesting would be a measurement of the x Pt spectrum, which is expected 
to be dominated by direct production rather than B decays. A preliminary measurement 
by CDF pH reports BR(ip -»■ /i+/i~)x cr(Xc -> ^7! Pt x > 7GeV; \rj\ < 0.5) = 3.2 ± 
0.3 ± 1.2 nb. Both xi an d X2 are here included. This can be compared with the range 
0.64nfe < a < b.lnb obtained using the LO QCD calculation described above ||119|| . Using 
the above cross section and using the inclusive B — > Xci branching ratio of 0.54 ± 0.21% 



128 ], we estimate that a fraction smaller than 10% of the %'s comes from B decays. 

A measurement of the production cross section and p t spectrum for T states would 
also be very useful in understanding the quarkonium production mechanisms. In this 
case one has at least three advantages: (i) the masses involved are larger and presumably 
both the non-relativistic approximation involved in the determination of the quarkonium 
wave function and QCD PT would work much more reliably than for charmonium; (ii) 
the signal does not have a contamination similar to the one due to B decays; (iii) the p t 
spectrum could be extended to very small values of p t , thanks to the large mass of the T 
and the large momentum of the decay muons. 



5 W AND Z PRODUCTION 

5.1 Inclusive Measurements 

Inclusive production of W and Z bosons is the most accurately known process in hadronic 
collisions. The absence of final state strong interactions affecting the observed state, one 
or two large-p t charged leptons, allows for high precision measurements and calculations. 
Uncertainties in the measurement of the total cross sections ||129| , |130|| are less than 10% 
and are dominated by the uncertainty on the absolute luminosity (see Table [TJ). 
The full NNLO 0(a s 2 ) corrections to the cross section are known ||131|| and techniques 



for the resummation of classes of leading and subleading logarithmic corrections to all 
orders of PT are available [ |132]| . The current theoretical systematic error is below 5%, 
estimated by varying factorization and renormalization scales within the range 10 GeV 
< fi < 1000 GeV. Slightly larger uncertainties arise from the use of different PDFs. The 
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Table 1: a w ■ BR, a z ■ BR and R = a w ■ BR/a z ■ BR at 630 and 1800 GeV. 
Data vs. 0(a s 2 ) QCD for different PDf sets pi] . BR(W -> ^)=0.109 and 
BR(Z -> l+l~)= 3.35xlQ- 2 . 





Data 


HMRSB 


MTE 


MTB 


a w • BR (pb) 630 GeV 


UA1 
UA2 


609 ±41 ±94 
682 ± 12 ± 40 


733 


699 


720 


ct z ■ BR (pb) 630 GeV 


UA1 
UA2 


58.6 ±7.8 ±8.4 
65.6 ±4.0 ±3.8 


69.2 


71.0 


69.9 


R (630 GeV) 


UA1 
UA2 


10.4l};| ±0.8 
10.4±g;g ± 0.3 


10.6 


9.9 


10.3 


(7 W • BR (nb) 1800 GeV 


CDF 


: 2.20 ±0.04 ±0.20 


2.06 


2.02 


2.10 


cj z • BR (pb) 1800 GeV 


CDF 


: 214 ± 11 ± 20 


194 


192 


198 


R (1800 GeV) 


CDF 


: 10.0 ±0.6 ±0.4 


10.6 


10.5 


10.6 



agreement between theory and experiment, at both SppS and Tevatron energies, is within 
one standard deviation and does not favour any particular set of PDF's provided one uses 
recent NLO fits. Even though the 0(a s 2 ) corrections add only a very small numerical 
contribution to the 0(a s ) result, they conspire to improve the stability of the cross section 



under changes of /x by a factor of 3-5, depending on the beam energy and PDF set [ 131 
This stability and the agreement with data represent a remarkable success of perturbative 
QCD. 

The charged lepton rapidity asymmetry in W decays: 

da/dy(£+) - da/dyjr) 
[V> da/dy{£+) + da/dy{£-)' 1 ' 

is more sensitive to the choice of PDF set and is not affected by luminosity uncertainties. 
Its measurement probes directly the quark components and the sea flavour symmetry of 
the proton ||133|| , necessary ingredients for a precise measurement of the W mass [134 



Current data at \y\ < 2 ||135|| already discriminate between different PDF fits. The 0(a s 



calculation of this asymmetry is available ||136|| , and new data will hopefully extend this 



measurement to more forward regions, where the difference between PDF's is expected 
to be even more pronounced. 

NLO calculations have also been recently completed for the inclusive W and Z p t 



distributions ||137|| . Measurements have been reported by UAl, UA2 and CDF ||138| , |139 
and are shown in Figure 120. The main source of systematic uncertainties in the case of the 
pY measurement is the determination of the neutrino transverse momentum, degraded 
by the energy resolution for the jets possibly present in the event. The small statistics 
(10% relative to the W case) limits instead the otherwise very clean pf measurement. 

The agreement with QCD is good at large p t , indicating consistency with the SM 
expectations. At smaller p t the theory is in better agreement with the UA2 data than 
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with CDF. The small pt region is interesting from the theoretical point of view, because a 
correct description of the spectrum requires the resummation of multiple gluon emission, 
which can be calculated in perturbative QCD | |140| | in the form of Sudakov form factors 
l2"fl . These effects have been included in the theoretical curves shown here ||141|| using 
the techniques developed in | |142| . Additional higher statistics measurements of the Z p t 
spectrum will help turning the qualitative agreement indicated here into solid QCD tests 
in the delicate semi-inclusive p t — > region. 



5.2 Associated Jet Production 

The production of jets associated with W's and Z's is less well predicted than the inclusive 
momentum spectra. Nonetheless, the characteristics of multijet final states in these events 
is very topical since it forms a background to top production. As with purely hadronic final 
states, most predictions for multijet characteristics in W and Z events are only available at 
tree level ||146|| , hence absolute cross section estimates have large uncertainties associated 
with the terms. Recent work has led to new NLO predictions for quantities such 
as the jet E t and pseudorapidity distributions in W+l jet events. 

CDF and UAl have measured the multiplicities of jets associated with W and/or Z 
production and have compared the results to tree level predictions ||143| , 144 1. Within the 
relatively large statistical and theoretical uncertainties the results are in good agreement 
with the theory. Figure [H| shows the cross section for W production as a function of jet 
multiplicity from CDF data. Other distributions, such as the E t distribution of associated 
jets, show some discrepancy with tree level predictions | |145| | . New NLO predictions will 



possibly improve the agreement with the data p4] . With more data at the Tevatron, it is 
expected that a more thorough test of W and Z plus jet production can be carried out. 

Using the ratio of the W+l jet and W+0 jet event rates, and comparing with the 
results of a LO calculation for W+l jet production, UA2 has extracted a measurement 
of a s (M^/): a s =0.123 ± 0.018( stat.) ± 0.017(syst.). This value is consistent with other 
determinations of a s from LEP and DIS data ||117|| . We point out that a fully consistent 
measurement of a s and an extraction of Aqcd can however only be performed using a 
NLO calculation for the W+l jet process. Only at this order it is possible to reduce the 
/i scale uncertainties and to define a precise renormalization scheme within which a s is 



measured. New analyses using the calculations of Ref. [24] will hopefully follow. 



6 DIRECT PHOTONS 

6.1 Single Photon Production 

As in the case of Drell Yan, the measurement of photons produced directly in a hadronic 
collision [|14S|| has the advantage of not suffering from final state strong interactions. Fur- 
thermore, since EM energy is detected with much better resolution than hadronic energy, 
systematic errors in the measurement of the photon momentum and direction are smaller 
than in jet measurements. Production of direct photons at small pt is dominated by pro- 
cesses with a qg pair in the initial state, be them of the Compton or of the bremsstrahlung 
type (Figure |22|). The capability of the experiments to observe direct photons at small p t 
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provides therefore yet another potential tool, in addition to the b quark measurements, 
to explore the gluon content of the proton at small values of x, or alternatively to learn 
more about small-x phenomena. The associate production of photons and charm quarks 
has also been suggested as a direct probe of the charm density in the proton ||149|| . 



Several difficulties however complicate the study of direct photons. First of all there 
are severe backgrounds to photon identification coming from hadrons such as ir° and 77's 
decaying into almost collinear photon pairs, faking a single 7. This background is statis- 
tically subtracted using two techniques. One technique relies on the different probability 
that one photon or a photon pair will convert in a e + e~ pair, the probability being in- 
dependent of pt- This "conversion method" can be used for arbitrarily large values of pt- 
The second technique relies on the measurement of the transverse shape of the EM shower 
in the calorimeter to determine the fraction of events with two overlapping photons. This 
"profile method" can only be applied over a limited p t range, above which the two photons 
are too close to be separable. 

On the theoretical side, predictions depend on the knowledge of bremsstrahlung con- 
tribution, which has both a perturbative and a non-perturbative piece. The latter is 
needed to properly define the boundary condition of the perturbative parton— > photon 
fragmentation function. It is due to the intrinsic hadronic component of the photon and 
it leads to a non-negligible g — > 7 fragmentation probability via Vector Meson Dominance 
(VMD). 

To reduce the hadron decay backgrounds, experiments do not measure a fully inclusive 
spectrum, but the so called isolated photon spectrum. Isolation is defined in different 
possible ways. UA2 requires no charged tracks within a A77 x A0 = 0.2 x 15° window 
around the 7 direction, and no EM energy within AR < 0.265. CDF requires the presence 
of less than 2 GeV of hadronic energy inside a cone of radius AR < 0.7 surrounding the 
photon. The isolation reduces the bremsstrahlung contribution ||153|| - fl55|l and emphasizes 
the purely perturbative effects, allowing for a more direct test of QCD. 

Full NLO calculations are available for the inclusive |[L52|| and isolated p t spectrum 



153, 157], as well as for the photon+jet processes [156]. A detailed study of the effects 



of isolation is presented in [ |154| |. The comparison between theory and data is shown in 
Figure |23|, which includes both UA2 and CDF results. While the agreement for p t > 20 
GeV is rather good, a discrepancy is apparent at smaller p t values. This is even more 
clear at the Fermilab energy. Several effects could be responsible for this problem. We 
will briefly survey them here. 

First of all, as always in PT, there is an intrinsic scale uncertainty. Here the scales 
needed are three: for renormalization, initial state factorization and final state fragmenta- 



tion. Studies reported in ||151|| indicate that the shape of the spectrum is rather insensitive 
to the scale uncertainty, at least in the p t range probed experimentally. Not even the use 
of different PDF sets can accommodate the factor of 2 discrepancy observed for the lowest 
p t bins ip.51]] . As in the case of the b cross section, Figure|T7"|, the values of p t are probably 
too large to allow significant departures from current PDF fits. 

The next possible effect is the bremsstrahlung contribution: how well do we know it? 
Ref. ][157|1 describes the full NLO correction to the bremsstrahlung processes, including 
a VMD description of the photon as a phenomenological input for the evaluation of the 
g — > 7 fragmentation. The results indicate that higher order terms add at most 50% to 
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the lowest order fragmentation contribution to the inclusive spectrum. After isolation 
cuts their effect will be even smaller, because the g — *> isolated-7 fragmentation is highly 
suppressed. We believe that 50% is therefore a reasonable estimate of the uncertainty 



reached today on the size of the bremsstrahlung contribution. Figure [24] shows the effect 
of removing the isolation requirement from the NLO QCD calculation [151 . This increases 
the QCD result by no more than 30%. A 50% uncertainty on this number is not sufficient 
to entirely explain the observed differences. 

We cannot exclude that a combination of all three effects just considered, in addition 
perhaps to new data and a better understanding of the experimental systematics, can 
reestablish agreement between theory and observations. Another possibility is however 
open. That is, the violation of naive factorization at small x, as was discussed in the heavy 
quark section. Like in that case, new diagrams with a t-channel gluon exchange appear 
at NLO for the first time (Figure ^2|). The same considerations and techniques outlined 
previously apply to this case | |109| ], even though no explicit calculation of the corrections 
to the differential pt spectrum has been carried out as yet. This issue will have to be 
properly understood before the photon distributions - either in p t or in rapidity - can 
be used to extract sensible measurements of the gluon structure functions in the small x 
region 



153, 158 



6.2 Double Photon Production 

Interesting measurements have also been performed on the direct production of photon 
pairs. Aside from its interest for QCD, this process is undergoing intense scrutiny as a 
possible dominant source of background to the detection of an intermediate mass Higgs 
boson at supercollider energies |]159||. The capability of QCD to properly estimate the 77 



production rate is therefore a very important fact to establish. 

Three processes contribute to the production of 7 pairs (Figure |25|): direct quark an- 
nihilation (qq — ■> 77, 0(a 2 )), gluon fusion via a quark box diagram (gg — > 77, 0(a 2 a 2 )) 
and various bremsstrahlung contributions (qg — > 577, 0(a 2 a s )). Even though of dif- 
ferent order in a s , these contributions are all comparable in magnitude over the cur- 
rently measured p t range, because at small x we have q(x) ~ a s g(x). The complete 
0(a 2 a s ) calculation is available [ |161|| , including the effect of isolation cuts |[L62j| , together 



with the LO gg — > 77 process. Data from UA2 ||150j| and CDF ||160|| are shown in Fig- 



ure |2B], compared to the relative calculations. In the case of UA2 the photons are not 
required to be isolated. Backgrounds and bremsstrahlung are reduced by applying the 
cut: p 4 (7i) ■ Ptilz) < ~ 0-7|Pi(7i)| 2 - The theory calculations reproduce the experimental 
selection criteria. 

The CDF data are systematically above the QCD curve, in particular at low p t . UA2 
shows a discrepancy only in the first p t bin. In addition to the pure QCD curve, the 
figure shows the results obtained by the PYTHIA shower MC, with and without the 
bremsstrahlung terms. The comparison between the different curves suggests that 1) 
PYTHIA has a bremsstrahlung contribution larger than NLO QCD and ii) initial state 
radiation induces a significant smearing of the p t spectrum. It is perhaps premature to 
formulate a judgement in relation to this measurement. On one side the statistical errors 
are still large. On the other the calculations have not been completed at the full 0(a 2 a s 2 ), 
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where we know some important contributions (gg — > 77) but we ignore the effect of others 
a priori comparable in size, such as gg — > qqjj- This last process would also contribute 
to a broadening of the 77 correlations w.r.t the available 0(a 2 a s ) estimates, which are 



unable to explain the data [|160|| . Last but not least, the values of x probed by this 



measurement are even smaller than those relevant for the b cross section, therefore this 
process is another interesting candidate for the study of small- 2 effects on production 
mechanisms. 

CDF also measures the average transverse momentum of the photon pair, (K t ) = 
5.1 ± 1.1 GeV. This is consistent with what expected from perturbative initial state 
radiation, (Kt) ~ a s (vl) ~ 4 GeV, considering that the bremsstrahlung processes will 
contribute an additional unbalance. CDF quotes agreement with the prediction of the 
PYTHIA calculation for (K t ). 



7 CONCLUSIONS 

All hard scattering processes in hadronic collisions require some understanding of QCD 
to be properly described. This is valid to the extent that they depend explicitly on a s 
and the parton distribution functions. Although QCD is widely accepted as the theory of 
strong interactions, progress can only result from successively making more rigorous tests, 
where discrepancies are not idley dismissed, but both data and theoretical assumptions 
are closely examined. 

In order to summarize the status of QCD predictions, one can imagine two ways of 
classifying results. In the first one could select phenomena according to the quality of 
the agreement between theory and experiment. In the second, one can select phenomena 
according to the presumed reliability of theoretical predictions and the corresponding 
faith in experimental results. It is a fact that processes which are believed to be reliably 
calculated also happen to belong to the class for which the agreement with data is good. 
This is the case for the 1-jet inclusive distributions and for W and Z production, which 
should be considered as successes of the application of perturbative QCD to hadronic 
collisions. There is, however, a possible discrepancy in xt scaling for jets, which should 
be an incisive test for the theory. As this article goes to print, there is no obvious 
explanation for such a discrepancy and we look forward to resolution, either via more 
data, or a new insight in the comparison to theory. 

In contrast to inclusive jet and W,Z production, there are processes such as b quark 
and direct photon production, where the theoretical uncertainties are large even at NLO. 
Perturbative K factors are big and strongly dependent on the choice of factorization and 
renormalization scales. Even worse, the disagreement between theory and data seems to 
be larger than the presumed uncertainties can account for. With independent data for 
parton distributions in this range of a:, it appears unlikely that one can find fault in a 
lack of knowledge of the gluon densities. There are, on the other hand, strong indications 
that a deeper understanding of the perturbative picture may be required to explain the 
discrepancies. In the case of b cross sections, more data, particularly with the power of 
secondary vertex detectors, will provide strong checks on the data. 

As indicated in the review of direct photon results, the processes contributing to 



23 



photon or heavy quark production at the NLO have singularities which are not present 
at tree level. For example, this is the case of diagrams with a t-channel exchange. Since 
these singularities only appear at NLO, an even higher order calculation would be needed 
to have a true NLO approximation to all relevant processes. This does not represent a 
problem for the 1-jet inclusive distributions or for the W, Z and Drell-Yan: in the first 
case no new singularity appears at NLO (t-channel gluon exchange is already there at tree 
level), in the second case, the available calculations are already at NNLO. This distinction 
could explain why there appears to be two classes of processes. 

Perturbative techniques for the study of multijet configurations are rapidly evolving 
and the agreement with data is quite reasonable. These tests are crucial to the search for 
new phenomena in events containing multiple jets. 

The measurement of finer details of the event structure, such as jet shapes, fragmen- 
tation and multijet correlations shows a good agreement with the results of both shower 
MC's and parton level calculations. This is therefore a success of perturbative QCD and 
of the way higher order processes are included in the MC algorithms. These measurements 
support the concept of local parton-hadron duality and establish a firmer ground for the 
use of shower MC's to predict the fine details of the jet structure in hadronic collisions. 
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FIGURE CAPTIONS 



Fig. 1: Gluon densities according to the most recent PDF analyses: MRSD0 [34] and 



various CTEQ fits [Pqj . For 0.01 < x < 0.1 and Q > 10 GeV, differences never 
exceed the 10% level. 

Fig. 2: Comparison of the inclusive jet cross section for LO QCD predictions with ex- 
periments at the ISR, SppS and Tevatron colliders |6f| |5l], [53|. Only one free 
parameter in the theory has been fixed (renormalization scale) in order to obtain 
this figure. 

Fig. 3: Examples of diagrams contributing to the jet cross section at C(af). Collinear 
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and soft singularities cancel between loop and tree diagrams, after imposition of 
a sensible jet definition involving finite opening angles for the final state partons. 



Fig. 4: 



Fig. 5: 



Fig. 6: 



Fig. 7: 



Fig. 
Fig. 



9: 



Sensitivity of the inclusive jet cross section to the renormalization scale for 0{a 2 s ) 
and predictions. Note that near the hard scattering scale, E t , the sen- 

sitivity is greatly reduced for the NLO calculation and the /i dependence goes 
from monotonic to forming a plateau near the hard scattering scale. 

The inclusive jet E t spectrum for CDF data using a cone size of 0.7, compared 
to theory as a ratio of (Data - Theory) /Theory. The upper plot (a) illustrates 
the theoretical uncertainty associated with variation of the renormalization scale 
fi (E t > fi > E t /4) for both LO and NLO. The lower plot (b) illustrates the 
dependence on the choice of PDF. The ) prediction using the HMRS set B 
PDF is used clS db reference. 



The variation of the jet cross section with clustering cone size R for jets of 100 
GeV E t . The standard 0(a^) calculation uses the merging parameter R sep = 2.0, 
whereas a modified version employs R sep = 1.3 fl67 . 



Fraction of energy contained in a sub-cone of radius r in jets found with a R=l 
cone. The data are from CDF charged tracking information, the QCD predictions 



are from |7] and HERWIG g|]. 



The ratio of dimensionless cross sections measured at \fs 
TeV compared to QCD predictions at both LO and NLO. 



1.8 TeV and 0.546 



The dijet angular distribution, dN/dx from CDF data shown along with 
0(a 2 s ) and 0(oil) predictions. The data are divided into three bins of dijet 
invariant mass, Mjj. 



Fig. 10: The dijet angular distribution from CDF data |76] for the highest values of Mjj 
compared with a model which includes both QCD and a parameterization for 
the effects of quark compositeness. 

Fig. 11: The evolution of the jet fragmentation function, F(z), as a function of dijet 
invariant mass, Mjj. This is shown along with fits of the form A InMjj + B |79| . 



Fig. 12: Angular separation for pairs of jets in four jet events. The solid line are the 
predictions from exact LO QCD matrix elements |]8*2"|| , and the dashed line repre- 
sents the expectations of phase space. The tree level predictions clearly describe 
the data much better than phase space. Jets are ordered by p t [j9Tj . 



Fig. 13: The distribution of jet pt for 4, 5 and 6 jet events from the UA2 collaboration 
|8"9"|| . The solid curve represents the exact LO QCD calculation |^2[ ^3] for four 
jets. The dashed-dotted line is the result of the Maxwell approximation [37 



for five jets, and the dashed lines are the predictions using the Kunszt-Stirling 
approximation 186 ]. 
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Fig. 14: The jet multiplicity plotted for different minimum jet pt cuts for events with 



greater than 400 GeV total transverse energy from CDF data | 93fl . The his- 
tograms are from the HERWIG event generator combined with a detector sim- 
ulation. Each histogram represents a different choice of PDF; employed were 
DOl (solid), D02 (short-dashed) 0, EHLQ1 (long dashed) and EHLQ2 (dot- 
dashed) pf. 



Fig. 15: A representative diagram for the t-channel gluon exchange contribution to heavy 
quark production. 

Fig. 16: Integrated b pt distribution at UA1 (left) and CDF (right): data versus NLO 
QCD. The lower curves correspond to (fi, Aqcd) = (mT,215 MeV), the upper 
ones to (/z, Aqcd) = (m T /4,275 MeV), with m T 2 = p t 2 + m Q 2 . 275 MeV cor- 
responds to one standard deviation from the central value of the MRSDO fit for 

A (2) 

Fig. 17: Fraction of the NLO QCD b cross section at 1.8 TeV coming from gluons with 
x g < x, for different p t thresholds. 

Fig. 18: A picture of the evolution of a gluon towards small-x. 

Fig. 19: J/ip Pt distribution at UA1 (left) and CDF (right): data versus QCD. Dotted 
line: direct quarkonium, dashed line: b decays, solid: total. The lower set of 
curves correspond to (/i, Aqcd) = 215MeV), the upper set to (/^Aqcd) = 



(%/ '4, 275 MeV). Parton distribution set MRSDO [3|| 



Fig. 20: W p t distribution at UA2 (right) and CDF (left): data versus QCD ililj . The 



band indicates the theoretical uncertainty due to the choice of factorization scale 
and PDF sets. 

Fig. 21: W+ n jet production rates at sfs = 1.8 TeV p3| . 

Fig. 22: Sample diagrams contributing to prompt photon production. Left: LO Compton 
scattering; Right: NLO bremsstrahlung. 

Fig. 23: Isolated prompt photon p t distribution at CDF and UA2, compared to a NLO 
QCD calculation ||153| |. For CDF, profile (circles) and conversion (diamonds) 
methods have separate normalization uncertainties, shown in the legend. 



Fig. 24: Study |151| of the effect of isolat ion on the photon p t spectrum at 1800 GeV. 



The solid lines indicate the relative variation of the theoretical calculation after 
reducing the isolation cone to 0.4, and after removing the isolation. 

Fig. 25: Sample diagrams contributing to double prompt photon production. 

Fig. 26: Double prompt photon p t distribution at 630 and 1800 GeV, compared to various 
theoretical calculations. The p t of both photons in each event enter in the plot. 
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